Study of the hadronization process in 
cold nuclear medium 

N. Akopo^|3, L. GrigoryaiJ^, Z. Akopov 
Yerevan Physics Institute, Br.Alikhanian 2, 375036 Yerevan, Armenia 

Abstract 

The improved two-scale model is used to perform the fit to the semi-inclusive 
deep-inelastic scattering (SIDIS) data of HERMES experiment at DESY on nuclear 
targets. The ratio of hadron multiplicity on nuclear target to the deuterium one is 
chosen as observable, as usually. The two-parameter's fit gives satisfactory agreement 
with the data in term of ^ criterium. Best values of parameters are then used 
to calculate the nuclear multiplicity ratio for the hadrons not included in the fit 
procedure. 
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1 Introduction 

Hadronic reactions in a nuclear medium, ei- 
ther cold or hot can shed additional light on 
the hadronization process. Numerous mea- 
surements of hadron production on nuclear 
targets in SIDIS of leptons [l]-[8] are avail- 
able. In ultra relativistic heavy-ion colli- 
sions the jet-quenching and parton energy- 
loss phenomena are observed [HI EHl- In 
each case the observed hadron yields are 
differed from those in the corresponding re- 
actions on free nucleons. In comparison 
with other reactions leptoproduction has 
the virtue that energy and momentum of 
the struck parton are well determined, as 
they are tagged by the scattered lepton. 
Study of hadron production in SIDIS on 
nuclear targets offers an opportunuty to in- 
vestigate the quark (string, color dipole) 
propagation in nuclear matter and the space- 
time evolution of the hadronization pro- 
cess. If the final hadron is formed inside 
the nucleus, it can interact via the rele- 
vant hadronic cross section, causing further 
reduction of the hadron yield. The per- 
turbative QCD cannot describe hadroniza- 
tion process because of the essential role of 
"soft" interactions. Therefore, the under- 
standing of this process on the phenomeno- 
logical level is of basic importance for de- 
velopment of the theory. For this purpose 
we investigate the nuclear attenuation (NA), 
which is a ratio of differential hadron mul- 
tiplicity on a nucleus to that on deuterium. 



where z = Eh/u, and u are energies 
of the final hadron and virtual photon re- 
spectively, N^{v,z) is the number of semi- 
inclusive hadrons at given v and z and N^{v) 
is the number of inclusive DIS leptons at 
given z/. Subscripts A [D) denote that re- 
action takes place on nucleus (deuterium) 
respectively. In the above formula more 



variables like the photon virtuality - Q^, 
transverse hadron momentum in respect to 
the virtual photon direction - pt , over which 
the NA is averaged, are not written. At 
present, several phenomenological models 
for description of the NA [TT]-|2B] are avail- 
able. The simple version of the string mo- 
del, so called Two-Scale Model (TSM), was 
proposed by European Muon Collaboration 
for the description of its experimental data |2]. 
In Ref . [I6j improved version of TSM (ITSM) 
was proposed. In present work ITSM is 
used to perform a fit to the recent SIDIS 
data of HERMES experiment on nuclear 
targets [7]. For a fit we use the more pre- 
cise (high statistic) part of the data sam- 
ple including one dimensional data for tt"*" 
and 7r~ mesons and two dimensional data 
for charged pions. One (two) dimensional 
data means that data are presented in form 
of function of one (two) variable. Then the 
R'I.j for all measured hadrons were calcu- 
lated with the values of parameters corre- 
sponding to the minimum values of reduced 

= X^/d-O-f. (here d.o.f. denotes "de- 
gree of freedom"). Then the results of such 
"best fit" were compared with the experi- 
mental data. 

The remainder of the paper is organized 
as follows. In section [2] we briefly remind 
about the ITSM. In section [3] the part of 
data included in fit and some details of fit- 
ting procedure are presented. Results are 
discussed and compared both with the dif- 
ferent versions of present fit and with our 
preceding one [TB|. In section H] we compare 
results of the fit with experimental data 
and discuss them. Conclusions are given 
in section [51 

2 ITSM 

Basic formula for NA in TSM [2j is: 

/■oo roo 

Ra = 2it bdb / dxp{b, x) x 

JO J-oo 

[1- / dx'a'''^{Ax)p{b,x')f-' , (1) 



where b is the impact parameter, x - longi- 
tudinal coordinate of the DIS point, p{b, x) 
- nuclear density function, x'- longitudinal 
coordinate of the string-nucleon interaction 
point, cr**^(Ax) - the string-nucleon cross 
section on distance Ax = x' — x from DIS 
point, A - atomic mass number. The above 
equation does not take into account the fi- 
nal state interactions (FSI) in deuterium. 
In this paper we use, following [16], more 
precise formula for the ratio of multiplic- 
ities R% = Ra/Rd- The string model is 
based on the idea that after DIS the knocked 
out (anti) quark does not leave the nucleon 
remnant, but forms a string (color dipole) 
with the (anti) quark on the fast and the 
nucleon remnant on slow ends, while the 
color string itself consists of gluons. Its 
longitudinal size must be larger than the 
transverse one, but cannot be essentially 
larger than the hadronic size because of 
confinement. The string can break down 
into two strings according to the following 
scenarios. First, when the quark-antiquark 
pair from the color field of the string is pro- 
duced; and second, when the color interac- 
tion between the string and the nucleon (ly- 
ing on its trajectory) has happened (see for 
instance [121 HH]). In the "history" of the 
string there are two time scales which are 
of interest to us. They are the time scales 
connected with the production of the first 
constituent (anti) quark of the final hadron 
and interaction of its two constituents for 
the first time. These two scales are (see 
Fig. H]): Tc (/c) - constituent formation time 
(length)!^]; and (Ih) - yo-yo formation time 
(length). The yo-yo formation means, that 
the colorless system with valence content 
and quantum numbers of the final hadron 
is formed, but without its "sea" partons. 
In the two-dimensional string model which 
satisfy the following conditions: (i) quark- 
antiquark pairs arising from the vacuum do 

''in relativistic units {h — c — 1, where h = 
h/2TT is the Plank reduced constant and c - speed 
of hght) Ti = li, i=c,h because partons and hadrons 
move with near hght speeds. 



X — DIS point 
t IL— const, formation length 

\,- yo-yo formation length 
L - full hadronizotion length 




X x+le x + l, X + L X 



Figure 1: Space-time structure of hadroni- 
zation in the string model. The two con- 
stituents of the hadron are produced at dif- 
ferent points. The first and second con- 
stituents of hadron h are created at the 
points Ps and P2, respectively. They meet 
at to form the hadron. 

not have energy; (ii) energy loss of the lea- 
ding quark on unit length (string tension) 
is constant (widely known example is the 
Lund model), there is a simple connection 
between and Tc 

Th-Tc = Zvj K , (2) 

where k - string tension (string constant). 
Further we will use two different expres- 
sions for Tc. The first expression is obtained 
for hadrons containing leading quark |29j : 

Tc = (1 - z)i'/k . (3) 

The color string fully spends its energy on 
the distance of L = v/ n beginning from the 
DIS point (see Fig.[T]). Last hadron produc- 
ing from the string is h = Hi, which con- 
tains leading quark and carries energy E^. 
At distance L, the energy of the leading 
quark becomes equal to zero and whole en- 
ergy of hadron is concentrated in another 
constituent. This constituent collects its 
energy from the string, and will have en- 



ergy Eh on distance L only if it was pro- 
duced on distance Eh/i^ = zu/k from L. 
This is reflected in eq. ([3]). It is important 
to note that the hadron produced on the 
fast end of string is not always necessarily 
the fastest hadron. Second expression for 
Tc used in this paper is its average value: 

r,= / ldlD,{L,zJ)/ / dlD,{L,z,l) , (4) 

where Dc{L, z, I) is the distribution of the 
constituent formation length I of summed 
over all ranks hadrons carrying momentum 
z. This distribution in framework of the 
standard Lund model [30] was obtained in 
Refs. [131 [31]: 

D,{L,z,l) = L{l + C)^^^^^x 
(.(/-L + .L) + i±|)x 

e{i)e{L- zL-i) , (5) 

where C = 0.3 is the parameter which con- 
trols the steepness of the standard Lund 
fragmentation function. The path traveled 
by the string between the DIS and interac- 
tion points is Ax = x' — x. In the TSM the 
string-nucleon cross section has form: 

a'*'{Ax)=9{T,-Ax)crg + 9{Th- 

Ax)e(Ax - T,)(Ts + 9{Ax - Th)ah , (6) 

where dg, Gs and an are the cross sections 
for interaction with the nucleon of the ini- 
tial string, open string (the string contain- 
ing first constituent (anti) quark of final had- 
ron on its slow end) and final hadron, re- 
spectively. In this model the string-nucleon 
cross section is a function which jumps in 
points Ax = Tc and r^. 

In reality the string-nucleon cross sec- 
tion starts to smoothly increase from the 
DIS point, and reaches the value of the 



hadron-nucleon one at Ax = r. 

Unfortunately, it is impossible to ob- 
tain cT'^*'" from perturbative QCD, at least in 
the region Ax ~ r. This means that some 
model for the shrinkage-expansion mecha- 
nism has to be introduced. In this work 
we use four versions of a^*''. Two of them 
having linear and quadratic dependence of 
the cross section on Ax/r, were taken from 
Ref . [32] . Let us briefly discuss the physical 
reason behind linear and quadratic depen- 
dence (see Ref. [33]). The QCD lattice cal- 
culations show that the confinement radius 
is much smaller than the mean hadronic 
radii. Consequently the color field in the 
hadrons is located in tubes with a trans- 
verse size much smaller than the longitudi- 
nal one. The valence quarks and diquarks 
are placed at the end-points of these tubes. 
In case of inelastic scattering, the interac- 
ting hadron-tubes intersect in the impact 
parameter plane. The probability of the 
crossing of the tubes is proportional to their 
length. This means that o"^*'" increases pro- 
portional to Ax/r. In the naive parton mo- 
del, the inelastic cross section of a hadron 
with a nucleon is proportional to the trans- 
verse area which is filled in by its partons, 
i.e. 0"^*'' increases proportional to (Ax/r)^. 
The first version of a**^ is based on quan- 
tum diffusion: 

a'*" (Ax) = e{T- Ax)[ag + {an - (Tg)x 

Ax/r] + e{Ax - r)ah , (7) 

where r = Tc + cAr, At = Th — Tc. We in- 
troduce the parameter c (0< c <1) in order 
to take into account a well known fact, that 
the string starts to interact with hadronic 
cross section soon after creation of the first 
constituent quark of the final hadron and 
before creation of second constituent. 

The second version follows from naive 
parton case: 

(J**" (Ax) = e{T - Ax)[ag+ 
(ah - (T,)(Ax/r)2] + e{Ax - r)a/, . (8) 



Two other expressions for o"^*'" were also 
used [II1[T1]: 

a^*^(Ax) = a, - (a, - a,)exp{ - ^) (9) 
and: 

a^'^iAx) = a, - (a, - a,)ea;j9( - (^)').(IO) 

One can easily note that at Ax/t <C 1 the 
expressions (9) and (10) turn into (7) and 
(8), respectively. At the first glance it may 
seems that the ITSM, as opposed to the 
TSM, is actually a 1-scale model. But one 
must note that r is a function of two scales 
T = {1 — c)tc + CTh whereas the parameter 
c regulates inclusion of each scale into r. 

3 Details of fit and 
results 

For the fit the semi-inclusive data [7] of 
HERMES experiment on four nuclear tar- 
gets (helium, neon, krypton, xenon) and 
deuterium were used. Only most precise 
(high statistic) part of data was used. It 
was consisting from two pieces: 

(i) the piece of the one dimensional data 
including the z/ - and z - dependences of 7r+ 
and 7i~ mesons. Each dependence consists 
from 9 experimental points, i.e. for this 
piece we have all together 144 points. The 
one dimensional data for nuclear multiplic- 
ity ratio are a functions of single variable 
u or z, whereas in model, R'l^ enters as a 
function of two variables u and z (the us- 
age of two variables allows one to avoid the 
problem of additional integration over z or 
u in eq.(l)). For this reason we introduce 
m the Rij, in case of one dimensional data, 
second variable by next way. In case of u 
- dependence, for each measured u bin the 
value of z (averaged over the given u bin), 
and in case of z - dependence, for each mea- 
sured z bin the value of j> (averaged over 
the given z bin) are taken from the exper- 
imental data; 



(ii) the piece of the two dimensional 
data, containing the charged pions data on 
the same nuclei. This part is available in 
form of detailed binning over u (z) and 
three slices over z (u). We would like to 
remind that slices over z are: first 0.2 < 
z < 0.4, second 0.4 < z < 0.7 and third 
z > 0.7 and over u: first 6 < u < 12GeV, 
second 12 < u < VJGeV and third 17 < 
V < 23.5GeV. Each slice of each depen- 
dence consists from 8 experimental points, 
besides third slices over u in z- dependence, 
which consist from 7 experimental points, 
i.e. in this piece we have 188 points. 

For one and two dimensional data we 
select all together 332 experimental points. 

Now let us turn to the discussion of 
the ingredients of the string model. One 
of the important parameters is the string 
tension (string constant) which determines 
the energy loss by leading quark on unit 
length. In this work it was fixed at a static 
value determined by the Regge trajectory 
slope [Ml [35] 

K = l/(27ray = IGeV/fm . (11) 

For calculations the following nuclear 
density functions (NDF) were used. For 
deuterium the hard core deuteron wave func- 
tions from Ref. [SB] were used. For ^ife the 
shell model [37] was used: 

A 2 (A - 4) r\ . r\ , , 
Pir) = Poij + -^^^)exp(-^),(12) 

where = 1.31 fm for "^He. For ^°iVe, 
^'^Kr and ^^^Xe the Woods-Saxon distri- 
bution was used 

p(r) = po/(l + exp((r - r^)/a)) . (13) 

Three sets of parameters for NDF's from 
eq.(13)were used for the fit: 
first set [3B], a = 0.54 fm , 
rA = (0.978 + 0.0206^1/3)^1/3 . (^^^-^ 

second set [39j, a = 0.54 fm , 

1 61 

TA = (1.19^1/3 -—^)fm; (15) 



a"*" (7) 


NDF 


cTq m6 


c 




1 


0.62 ±0.15 


0.247 ±0.017 


0.63 


2 


0.81 ±0.16 


0.231 ±0.019 


0.57 


3 


0.91 ±0.16 


0.202 ±0.017 


0.61 


(t"*^'(8) 


NDF 


(Tq mb 


c 




1 


3.71 ±0.10 


0.175 ±0.010 


0.55 


2 


3.90 ±0.10 


0.161 ±0.009 


0.50 


3 


4.12 ±0.12 


0.140 ±0.014 


0.57 


a"*'' (9) 


NDF 


CTq mb 


c 


x' 


1 


1.29 ±0.12 


0.079 ±0.012 


0.55 


2 


1.53 ±0.13 


0.065 ±0.012 


0.53 


3 


1.76 ±0.14 


0.040 ±0.012 


0.59 


a"*" (10) 


NDF 


(Tq mb 


c 


x' 


1 


3.91 ±0.10 


0.100 ±0.012 


0.56 


2 


4.13 ±0.11 


0.087 ±0.012 


0.53 


3 


4.39 ±0.12 


0.066 ±0.012 


0.60 



Table 1: Values of fitting parameters and 
in case of Tc{3) and total errors. For 
NDF from eq.(13) versions 1, 2, 3 are sets 
of parameters from eqs.(14), (15), (16), re- 
spectively. a^^'^{7) means o"'**'" from eq.(7) 
etc. 

third set [30], a = 0.545 fm , 

ta = 1.14A^/^ /m . (16) 

The corresponding values of po were deter- 
mined from the normahzation condition: 

j d^rpij-) = 1 . (17) 

Parameter a is practically the same for all 
three sets, radius ta for the third set is 
larger by approximately 6 % than the ones 
for the first and second sets. Let us briefly 
discuss the choice of the nuclear matter dis- 
tribution functions. For deuterium the choice 
of the NDF is not important because the 
FSI are small. For light nucleus '^He the 
shell model was used, because there was 
no alternative. For middle and heavy nu- 
clei preferable NDF is Woods-Saxon distri- 
bution. However, there is some freedom 



a''' (7) 


NDF 


Gq mb 


c 


x' 


1 


0.00 ±0.01 


0.372 ±0.012 


1.89 


2 


0.00 ±0.01 


0.343 ±0.011 


1.71 


3 


0.00 ±0.01 


0.313 ±0.012 


1.67 




NDF 


cTq mb 


c 




1 


1.27 ±0.11 


0.146 ±0.004 


1.01 


2 


1.41 ±0.22 


0.137 ±0.018 


0.91 


3 


1.28 ±0.23 


0.113 ±0.016 


0.92 


CT^*^(9) 


NDF 


Gq mb 


c 


x' 


1 


0.00 ±0.01 


0.111 ±0.009 


0.97 


2 


0.00 ±0.01 


0.085 ±0.009 


0.85 


3 


0.00 ±0.01 


0.057 ±0.009 


0.84 


a^*'"(10) 


NDF 


(jq mb 


c 


x' 


1 


1.82 ±0.14 


0.077 ±0.010 


0.74 


2 


1.99 ±0.14 


0.068 ±0.010 


0.67 


3 


2.03 ±0.16 


0.047 ±0.010 


0.71 



Table 2: Values of fitting parameters and 
in case o/rc(4) and total errors. 

in the choice of the parameters themselves, 
therefore we have included three sets of pa- 
rameters (14)-(16), in order to study uncer- 
tainty of the fitting procedure related to the 
NDFs. Two expressions for were used for 
the fit - equations (3) and (4). For a^^^{Ax) 
four different expressions from eqs.(7)-(10) 
were used. The values of ah (hadron-nucleon 
inelastic cross section) used in the fit were 
set equal to: 0"^+ = 0",^- = 20 mb. The 
same value of inelastic cross section was 
used for charged pions. 

The fit was performed to tune two para- 
meters: the initial value of string-nucleon 
cross section cxg and coefficient c. 

The quantitative criterium was used. 
As usually it was determined as: 



a"*" (7) 



NDF 


CTq mb 


c 




1 


0.87 ±0.04 


0.256 ±0.007 


9.34 


2 


1.07 ±0.03 


0.241 ±0.006 


7.47 


3 


1.19 ±0.04 


0.202 ±0.006 


9.32 



NDF 




c 




1 


3.87 ±0.02 


0.165 ±0.004 


8.49 


2 


4.07 ±0.03 


0.150 ±0.005 


6.55 


3 


4.32 ±0.03 


0.123 ±0.005 


9.06 



a"*'' (9) 



NDF 


cTq m6 


c 


x' 


1 


1.49 ±0.03 


0.066 ±0.004 


8.31 


2 


1.73 ±0.03 


0.052 ±0.004 


7.02 


3 


1.99 ±0.03 


0.020 ±0.004 


9.15 



a"*" (10) 



NDF \ (Tq mb \ c \ 

1 4.05 ±0.02 "01)81 ±0.004 8.91 

2 4.26 ±0.02 "01)68 ±0.004 7.36 
~^ 4.56 ±0.03 0.042 ±0.004 9.80 



CT^*''(7) 



NDF 




c 




1 


0.00 ±0.00 


0.526 ±0.005 


38.3 


2 


0.00 ±0.00 


0.486 ± 0.004 


34.8 


3 


0.00 ±0.00 


0.441 ± 0.004 


34.8 



o-**'"(8) 



NDF 


cTq m6 


c 




1 


0.93 ±0.04 


0.141 ±0.004 


21.6 


2 


1.14 ±0.04 


0.141 ±0.004 


19.6 


3 


0.83 ±0.05 


0.099 ± 0.004 


20.2 



(T^*'(9) 



NDF 




c 


x' 


1 


0.00 ±0.00 


0.157 ±0.003 


21.4 


2 


0.00 ±0.00 


0.121 ±0.003 


18.2 


3 


0.00 ±0.00 


0.079 ± 0.003 


18.5 



a^*'"(10) 



NDF \ aq mb \ c \ 

1 1.54 ±0.03 0.059 ±0.003 15.8 

2 1.73 ±0.04 0.053 ±0.003 13.9 

3 1.70 ±0.04 0.023 ±0.004 15.3 



Table 3: Values of fitting parameters and 
X^ in case of Tc{3) and statistical errors. 

ARUexp) ) ' 

where riexp and ripar are numbers of exper- 
imental points and parameters; R\.[{theor) 
is the theoretical value for ratio at given 
point; R\.[{exp) and Ai?^(exp) are exper- 
imental values of R\,i and its error. Let us 
firstly discuss fit with use of total errors. 
Results are presented in the Tables [Hand El 
Easily to see that ITSM describes data on 
quantitative level. The version with Tc{ [3]) 
gives for x^ the values in order of 0.5 for 
all versions of a'**'". The version with Tc{^ 
gives for x^ the values in order of unity. 
The common tendency is that a**'' is es- 
sentially smaller than hadron-nucleon cross 
sections and coefficient c essentially smal- 
ler than unity. The obtained values of x^ 
show, that total errors are still large and 
do not allow to verify the different versions 
of ITSM. It is a reason why we turn to 
the fit of the data with statistical errors 



Table 4: Values of fitting parameters and 
X^ in case of T(.{^) and statistical errors. 

only. Corresponding results are presented 
in the Tables [3] and HI Again, as in case 
of total errors the values of x^ for Tc{ [3]) 
approximately two times smaller than for 
Tc ( HD but now values of x^ larger as min- 
imum ten times in comparison with pre- 
ceding case. For version with in form 
of eq.(3) minimum x^ = 6.55 is obtained 
in case of a**'' corresponding to eq.(8) and 
NDF in form of eq.(15). For version with 
Tc in form of eq.(4) minimum x^ = 13.9 
is obtained in case of a**'" corresponding to 
eq.(lO) and NDF in form of eq.(15). At the 
first glance x^ are too large and it is impos- 
sible to speak about quantitative descrip- 
tion of data. But we would like to remem- 
ber that statistical errors of data obtained 
by HERMES experiment and using for this 
fit [7j are in order of 1 % and for some 
points even 0.5 %. For experimental points 
measured with such small errors we obtain 
X^ = 9 when difference between theoreti- 



cal model and experimental value is 1.5 % 
only. We think, that = 6.55 which is ob- 
tained in case of rc(3), a'^''{8) and NDF(15) 
is suitable for description of experimental 
data. The agreement of version with Tc(4) 
is worse. We will discuss comparison of 
these two versions in the next section. 

Except the basic fit which was discussed 
above, three another fits were also performed. 
First of them was performed with the same 
data set as basic one but three parameters 
were used. Two of them were taken the 
same as in basic fit and as a third para- 
meter the string tension k, was used. The 
results of fit with Tc{^ were obtained very- 
close to the basic fit for both fitting pa- 
rameters and x^. Indeed, for (j**''(8) and 
NDF in form of eq.(15) the values cxq = 
3.87±0.042m6, c = 0.18±0.01, k = 1.075± 
O.OlSGeV/ fm and = 6.54 were obtained. 
In this case inclusion of third parameter 
does not felt. The results of fit with Tc{ H]) 
show the stronger dependence from the third 
parameter. For this case the following val- 
ues for cr**^(10) and NDF(15) were obtained: 
aq = 3.36 ± 0.058m6, c = 0.00 ± 0.00, 
K = 0.73±0m8GeV/fm and = 11.7. It 
is easily to see that the values of parameters 
significantly differ from the corresponding 
values in basic fit, but x^ is close enough 
to the one for basic fit. Two other fits were 
performed separately for one and two di- 
mensional pieces of data. Results are close 
to the ones for basic fit. When data are di- 
vided on one and two dimensional parts, 
the x^ for one dimensional data are ob- 
tained slightly larger than for two dimen- 
sional part because in first case statistical 
errors are smaller. But general situation is 
very resemble to the one for basic fit. 

Let us compare the procedure and re- 
sults of present fit with the our preceding 
fit [16]. At that time the SIDIS data of 
HERMES experiment for two targets: ni- 
trogen |4j and krypton [5] were available. 
The z/- and ^-dependences of vr"*" and vr" 
mesons with 58 available experimental points 
were included in the fit procedure. As the 



experimental errors only statistical ones were 
taken. Minimum values of x^ (best fit) for 
ITSM two parameter's fit were obtained for 
Tc from eq.(3) {x^ = 1.4) and for from 
6q-(4) (x^ = 1.5). Comparison with the 
new fit showes that x^ became several times 
larger. It should be pointed out that rea- 
sons for this increase are: (i) very small val- 
ues of statistical errors of present data (as 
a minimum two times smaller than of pre- 
ceding data); (ii) in the preceding fit were 
included data for two nuclei only (nitro- 
gen and krypton). The data on nitrogen 
had smaller number of experimental points 
and essentially larger errors than data on 
krypton, i.e. fit was mainly based on kryp- 
ton data. In present case data are available 
for four nuclei from light to heavy (helium, 
neon, krypton and xenon) and these data 
have comparable statistics. 

4 Comparison with data 
and discussion 

The results of the performed fit are pre- 
sented in Tables [H El |3] and |2 Tables [U 
and |2] show the best values of the fitted pa- 
rameters, their errors and x^ for the case of 
total errors. Tables [3] and S] show these val- 
ues for the case when only statistical errors 
were used for the fit procedure. Further 
we will discuss only last case. In Figs. |2] 
and [3] the one dimensional data for the 
hadron multiplicity ratio R^.f for vr^ and 
7r~ mesons as a function of variable u (left 
panels) and z (right panels) are presented. 
The results for ^ife(panels a and b), "^^Ne 
(c, d), ^^Kr (e, f) and ^^^Xe (g, h) are 
presented. Experimental points were taken 
from Ref. [7] (filled circles). For compari- 
son we presented in case of krypton previ- 
ous HERMES data [5] (open circles) also. 
The theoretical curves were calculated with 
the best values of parameters obtained for 
the versions of model with constituent for- 
mation length in form of eq.(3) (dashed 
curves) and eq.(4) (solid curves). Version 
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Figure 2: One dimensional data. Hadron 
multiplicity ratio R\j for tt"*" mesons as a 
function of variable v (left panels) and z 
(right panels). The results are presented 
for ^He (a, h), ^^Ne (c, d), ^^Kr (e, f) 
and^^^Xe (g, h). Experimental points from 
Ref. ^ (filled circles) and Ref. (open 
circles). The curves were calculated with 
the best values of parameters obtained for 
the constituent formation length in form 
of eq.(3) (dashed curves) and eq.(4) (solid 
curves). These data were included in fit. 

with Tc in form of eq.(3) describes data pre- 
sented in Figs. [2] and [3] better than ver- 
sion with Tc in form of eq.(4), but size of 
presented figures does not allow to see this 
difference. 

In Fig. m the two dimensional data for 
R'li for charged pions on '^He (panels a, b) 
and "^^Ne (panels c, d) functions of u 
(left panels) and z (right panels) are pre- 
sented. In Fig. Elthe two dimensional data 
for R'Ij for charged pions on ^'^Kr (panels 
a, b) and ^^^Xe (panels c, d) as a functions 
of u (left panels) and z (right panels) are 
presented. Experimental points are taken 
from Ref. [7]. Filled symbols are: triangles 
- experimental points for the first slice over 
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Figure 3: The same as described in the 
caption of the Fig. 2 done for n~ mesons. 
These data were included in fit. 

z (u) ioT u - {z - ) dependence; circles - for 
the second slice and stars - for the third 
slice. The open symbols represent the one 
dimensional data and were included in fig- 
ures for comparison with two dimensional 
ones. The open circles were chosen for 7r+ 
and open triangles for tt~ mesons. Unfor- 
tunately the two dimensional data for he- 
lium are not informative because the points 
from different slices are mixed. Data for 
"^^Ne are more useful because the separa- 
tion of points from different slices partly 
takes place. The complete separation of 
data from different slices takes place for 
heavy nuclei (krypton and xenon). They 
carry additional information in comparison 
with one dimensional data and are impor- 
tant for the development of theoretical mod- 
els. From these figures we see, that one 
dimensional data (open points) mainly co- 
incide with second slices and do not reflect 
the behavior of data in first and third ones. 
The theoretical curves were calculated with 
the best values of parameters obtained for 
the constituent formation length in form 
of eq.(3) (dashed curves) and eq.(4) (solid 




Figure 4: The two dimensional data. The 
ratio R'Ij for charged pions on ^He (panels 
a, h) and ^°iVe (c, d) nuclei as a function 
of u (left panels) and z (right panels) . Ex- 
perimental points from Ref. Filled sym- 
bols are: triangles - for the first slice; cir- 
cles - for second slice and stars - for third 
slice. Open symbols represent one dimen- 
sional data: circles - for 7r+ and triangles - 
for TT^ mesons. The theoretical curves were 
calculated with the best values of parame- 
ters obtained for the constituent formation 
length in form of eq.(3) (dashed curves) 
and eq.(4) (solid curves). These data were 
included in fit. 

curves). Let us remind that Tc of eq.(3) 
was obtained for leading hadron while Tc of 
eq.(4) is the average value of this quantity 
in the standard Lund model. As we already 
mentioned above the two dimensional data 
for helium presented in Fig. H] are less use- 
ful for comparison of different versions of 
model. We can state only, that theoret- 
ical curves do not contradict data. For 
middle nucleus (neon) presented in Fig. HI 
despite on partial mixing of experimental 
points from different slices, we can state, 
that version with Tc corresponding to eq.(3) 
describes data better than one with Tc cor- 



Figure 5: The same as described in the cap- 
tion of the Fig. [7] done for ^^Kr (panels a, 
b) and ^^^Xe (c,d) targets. These data were 
included in fit. 

responding to eq.(4). And, at last, let us 
turn to the heavy nuclei (krypton and xenon) 
where situation is more clear. From Fig. |5] 
we see that in average the version with Tc 
for leading hadron describes two dimensional 
data better than one from standard Lund 
model. The leading hadron approach satis- 
factory describes all three slices over z for 
p - dependence (excepting last points), and 
also the first and second slices over u for z 
- dependence but underestimates data in 
third slice, although the behavior is true. 
The main problem of the version of model 
with Tc taken from standard Lund model 
(eq.(4)) is the unsatisfactory description of 
first and second slices over z in u - depen- 
dence. This version of model practically 
does not differ data with small and middle 
z. The z -dependence it describes on satis- 
factory level excepting the region of small 
z. It is worth to mention that the difficul- 
ties in the description of z/ - and z - depen- 
dences are mutually connected. 

Furthermore, the NA for one dimen- 
sional data of hadrons produced on all nu- 




Figure 6: One dimensional data. Had- 
ron multiplicity ratio R\.^ for ii^ mesons 
as a function of variable v (left panels) 
and z (right panels). The results are pre- 
sented for ^He (a, b), ^^Ne (c, d), ^^Kr (e, 
f) and "^^^Xe (g, h). Experimental points 
were taken from Ref. [7J (filled circles) and 
Ref. [5^] (open circles). The curves were 
calculated with the best values of parame- 
ters obtained for the constituent formation 
length in form of eq.(3) (dashed curves) 
and eq.(4) (solid curves). These data not 
included in fit. 

clear targets (but not included in fit), were 
calculated in our model. In Figures 6-9 
the u - and z - dependences are presented 
for 7r°, K~^, K~ mesons and antiprotons, 
respectively. The hadron multiplicity ratio 
R\j for mentioned hadrons as a function of 
variable v (left panels) and z (right panels) 
are presented for ^He (a, b), '^^Ne (c, d), 
^"^Kr (e, f) and ^^^Xe (g, h). Experimen- 
tal points were taken from Ref. (filled 
circles). For ^'^Kr were taken also data 
from Ref. [5] (open circles). The theoretical 
curves were calculated with the best val- 
ues of parameters obtained for the versions 
with Tc in form of eq.(3) (dashed curves) 



Figure 7: The same as described in the cap- 
tion of the Fig. done for mesons. 
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Figure 8: The same as described in the cap- 
tion of the Fig. done for K~ mesons. 

and eq.(4) (solid curves). The following 
values of inelastic cross sections ah were 
used: cr^o = ax- = 20 mb, aK+ = 14 mb 
and (jp = 42 mb. Calculations, which were 
performed without additional fit, satisfac- 
tory describe data for vr^ (Fig. [6]) and K~ 
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Figure 9: The same as described in the cap- 
tion of the Fig. done for antiproton. 

(Fig. IH]) mesons. Here we also have, as 
in case of vr"*" and vr" mesons, that ver- 
sion with Tc in form of eq. (3) describes data 
shghtly better than with Tc in form of eq. (4) . 

It is worth to discuss data for mesons 
from Fig. [7] more in details, because here 
we have some problems. Let us begin with 
u - dependence. The data on ^He are de- 
scribed very well. For "^^Ne nucleus we have 
disagreement between experimental data and 
theoretical model for the first two points. 
And, at last, for heavy nuclei we have es- 
sential underestimation of data by theoreti- 
cal model. In this case comparison with the 
previous data for krypton helps to under- 
stand that problem rather in model than 
in data. In case of z - dependence we have 
satisfactory agreement for '^He and "^^Ne. 
Again description is worse for heavy nu- 
clei. First three points are underestimated 
by model. We see that in region of small z 
experimental values of R% decrease, while 
theoretical ones are rather constant. There 
is some disagreement between experimen- 
tal data and theoretical model in case of 
antiprotons also (see Fig. |9]). In this case 
theoretical curves slightly overestimate data. 



But we do not discuss it here in detail be- 
cause of two reasons: (i) large errors and 
(ii) some differences of the antiproton pro- 
duction mechanism from mesons one, which 
do not taken into account in present paper. 

5 Conclusions. 

The recent HERMES data [7J were used to 
perform the fit for ITSM. Two-parameter's 
fit demonstrates satisfactory agreement with 
data. The main goal of this paper was the 
further development of the model, in par- 
ticular the choice of proper version for r^. 
Minimum (best fit) was obtained for ver- 
sion of model with Tc in form of eq.(3). The 
version with Tc in form of eq.(4) gave es- 
sentially larger value for x^. Comparison 
with the two dimensional data obtained, 
for the first time, by HERMES experiment, 
allowed to perform additional verification 
of the different choices of constituent for- 
mation length. In particular it was ob- 
tained, that version with in form of eq. (3) 
has difficulty in description of third slice 
over u in z -dependence, i.e. in region of 
large u. The difficulties of version with Tc 
in form of eq.(4) are more serious, because 
this version does not differ small and mid- 
dle z in u - dependence. Although version 
of model with in form of eq. (3) describes 
data on the satisfactory level, the two di- 
mensional data show, that we have some 
problem connected with the choice of the 
more adequate form for constituent forma- 
tion length. More detail two dimensional 
data for identified hadrons that is expected 
from HERMES will provide essentially bet- 
ter conditions for the choice of preferable 
version of the model in terms of different 
expressions for and, may be, o"^*^. In all 
versions we have obtained that aq <^ ah- 
This indicates that at early stage of had- 
ronization process the color transparency 
takes place. It is worth to mention that 
the same result was obtained in our pre- 
ceding fit [TB]. We do not include in con- 



sideration the NA of protons, because in 
this case additional mechanisms connected 
with color interaction (string-flip) and final 
hadron re-scattering become essential (see 
for instance Refs. [121 113]) • 
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